New Approach to Continuum Path Integrals for Particles and Fields 
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An approach to approximate evaluation of the continuum Feynman path integrals is developed 
for the study of quantum fluctuations of particles and fields in Euclidean time-space. The paths are 
described by sum of Gauss functions and are weighted with exp(— 5") by the Metropolis method. 
The weighted smooth paths reproduce properties of the ground state of the harmonic oscillator in 
one dimension with more than about 90% accuracy, and the accuracy gets higher by using smaller 
width of the Gauss functions. Our approach is applied to quantum field theories and quantum 
fluctuations of U(l) and SU(2) gauge fields in four dimensions respectively provide the Coulomb 
force and confining linear potential at qualitative levels via the Wilson loops. Distributions of large 
values of gauge fields are found to be suppressed at least exponentially. 
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It is quantum physics which dominates microscopic 
phenomena less than the atomic scale [l|. One elegant 
way to describe the quantum phenomena is the path in- 
tegral method developed by Feynman Q, in which all 
possible paths are taken into account with the probabil- 
ity amplitude exp(iS/h) with S the action of the system. 
The path integral method gives a clear interpretation of 
behavior of particles in quantum mechanics as fluctua- 
tions from the classical paths, although exact evaluations 
of the path integrals are possible only in few cases such 
as the harmonic oscillator 0] . The path integral method 
also promotes modern developments of theoretical ele- 
mentary particle physics and supplies a nonpcrturbativc 
technique for quantum field theories Q. 

The evaluation of the path integrals can be simplified 
by discretizing time-space, in which derivations and inte- 
grations are replaced with finite differences and summa- 
tions, respectively, and measure of the path integrals be- 
comes a countable product. An important example of the 
discretized path integrals is the lattice QCD (quantum 
chromodynamics) [f|, by which nonperturbative aspects 
of QCD have been revealed @. However, the time-space 
discretization explicitly breaks continuous symmetries of 
time-space such as the translational symmetry down to 
discrete symmetries, and sometimes leads to qualitative 
discrepancies such as magnetic monopoles in the lattice 
QED (quantum electrodynamics) Q- Hence it is desired 
to perform the path integrals in continuous time-space 
from viewpoint of complementarity for the discretized 
approach. 

In this paper we develop an approach to evaluating the 
continuum path integrals in Euclidean time (i — > —it) for 
particles and fields. The continuous paths are described 
by sum of smooth functions with weight cxp(— S) by the 
Metropolis method Here we take the natural units 



H = C = = 1. 

Firstly, for a nonrclativistic particle with one degree 
of freedom q in a periodic boundary condition with pe- 
riod T, q(r + T) = q(T~), the path integral method evalu- 
ates the quantum transition amplitude in Euclidean time 
as 0, 



Z = 



period 



Vqexp(-S[q}), Vq = \{dq(T). (1) 



Here the measure T>q is formally defined as an uncount- 
able product, and the expression ([T]) means that the 
quantum transition amplitude corresponds to the sum- 
mation of all possible paths for the particle with the 
probability amplitude exp(— S). Since the quantum fluc- 
tuations of the particle are weighted with the factor 
cxp(— S), an expectation value of an operator 0[q] in 
quantum mechanics can be evaluated by using N paths 
q n (n = 1, 2, • • • , N) weighted with exp(— S) as, 
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VqO[q] exp(-S) 
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71 = 1 



where the last approximation becomes good for large N . 

Before explaining our approach to the continuum path 
integrals, we briefly review the discretized approach to 
the path integrals for a nonrclativistic particle in the pe- 
riodic boundary condition. In this case the particle posi- 
tion is represented as qi at time fj = ia with i = 1, • • • , 
Mat and the lattice spacing a = T /N\ atl and the path of 
the particle is obtained by connecting qi and qi+\ from 
i = 1 to iVi a t with straight lines. Then measure of the 
path integral is defined as a countable product, 



Vq = \\ dqi 



(3) 
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and the action is replaced with the corresponding dis- 
cretized one, S[q\. The discretized path integrals are 
evaluated in the following way (see 0]). Namely, change 
of the particle position at each time f,-, dqi, is generated 
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as a random number within [—A, A] with a fixed value 
A, in which Sqi has a property of the micro- reversibility. 
This change Sq, is judged by the Metropolis test [1], in 
which qi + Sqi is redefined as qi in acceptance probability 
min[l, exp(S[q]~ S[q + Sq])} and otherwise Sqi is rejected. 
We denote this step as Wi. Then the whole positions are 
updated by the "sweep", i.e., performing W% from i = 1 
to Mat- After several sweeps quantum paths in equilib- 
rium weighted with exp(— S) are obtained. 

Now let us make our procedure for the simulation of the 
continuum path integrals from analogy to the discretized 
approach. Our idea here is to connect the neighboring 
points for the particle positions by smooth lines rather 
than straight lines. In order to achieve this, we smear 
the particle position in the discretized approach qi with 
a Gauss function of width £j at time r, , 



q~i (at r = n) 



qi exp 



(4) 



for i — l to Mat, where r — t, means to take time dis- 
tance between r and tj in the periodic boundary condi- 
tion. Here (q^ Tj, in the continuous approach corre- 
sponds to ((/j, fj, a) in the discretized approach. With 
this smearing one can naturally connect the path with 
smooth lines rather than straight lines as, 



£ 

t=i 



q% cxp 



(5) 



where N sum (= Mat) is number of the summed terms. In 
this construction of the smooth path, as qi for each i takes 
value in range [—A, A] with a some cut-off A according 
to Eq. ([3]), q(r) at each time takes value in similar range 
[—A', A'] with a cut-off A' similar to A, and the minimal 
scale of the fluctuation for q(r) corresponds to The 
positions of the Gauss functions r, may take random val- 
ues rather than values in same interval, tj = iT/N sum , as 
long as every time is equally treated without making any 
special time. The width of the Gauss function, or the 
scaling constant, ^, is fixed so that £j does not depend 
on n, which prevents any special time. In this study 
we take two strategies for one is to generate £, ran- 
domly within [Aj, A^] in uniform probability with ultra- 
violet and infrared cut-offs A^ and A^, respectively (ran- 
dom scale) , and the other is to use a properly fixed value 
(fixed scale). Then the weight exp(— S) is given by the 
Metropolis method for qi as in the discretized approach. 
Namely, an additional fluctuation Sqi with randomly cho- 
sen i is randomly generated within [— A 9 , A q ] in uniform 
probability, where A q is a cut-off for the fluctuation am- 
plitude. Then the charge 5qi with respect to qi is judged 
by the Metropolis test, in which qi + Sqi is redefined as qj 
in acceptance probability min[l, exp(S[q] — S[q+Sq])] and 
otherwise 5qi is rejected. We emphasize that the addi- 
tional fluctuation Sqi is micro-reversible without making 
any special directions in coordinate space. 

Here we note that q(r) in Eq. ([5]) cannot describe all 
possible paths, since the Gauss functions in Eq. © can- 



not be a complete set with respect to the smooth func- 
tions in the periodic boundary condition. Therefore, at 
this point our construction of paths is an approximation 
with respect to the complete paths required by the mea- 
sure Vq. Nevertheless, we expect that most of the possi- 
ble paths can be taken into account when value of iV sum 
is sufficiently large with large A g , small A^, large Af , and 
dense n, in which case distribution of the Gauss functions 
with various width is sufficiently dense. We also note that 
time uniformity may be broken when the time compo- 
nents Ti take values in same interval, Tj = iT/N aum , but 
we expect that the uniformity will restore if one considers 
sufficiently dense distribution of the Gauss functions. 
Our procedure can be summarized as follows: 

1. Make an initial condition for the path by de- 
termining constants Tj, £j) from i = 1 to N snm 
in the following manner. Namely, n is generated 
within range [0, T] so as not to make any special 
time, and £i is randomly generated within [A^, A^] 
in uniform probability or is fixed as a proper value. 
The coefficient qi is randomly generated within 
[— A g , Aq] (hot start) or is taken to be zero for all 
i (cold start). 

2. Randomly choose i and generate Sqt within 
[— Aq, Aq] in uniform probability so as to construct 



an additional fluctuation, 



Sq(r) = Sqi exp 



(r n) 



2 1 



(6) 



3. According to the Metropolis method 0], ac- 
cept the additional fluctuation 8q in probability 
min[I, cxp(S , [q , J — S[q + Sq])]. If and only if the 
additional fluctuation Sq is accepted, we redefine 
the path q + Sq as q (or cquivalently redefine the 
coefficient qi + Sqi as qi). 

4. Iterate steps 2. and 3. until the action as well as 
other expectation values converge. 

In this procedure, we eventually obtain a smooth path for 
the particle ([5]), which is weighted with exp(— S) due to 
the step 3. Generalization to the systems with / degrees 
of freedom (/ > 2) is obvious. Here we emphasize that 
the description of smooth fluctuations only by the finite 
terms of Gauss functions is an approximation. 

Now let us examine our approach by investigating a 
harmonic oscillator in one dimension, which action is 
written as, 
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(7) 

with q = dqjdr. Here we fix its mass and angular fre- 
quency as in = uj = I, and we take the time range T = 
20. We consider three simulation conditions; A sum = 50, 
Aq = 3, and the fixed scale with £ = 1 (A), N sum = 100, 
Aq = 3, and the random scale with = 0.2 and A^ = 1 
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FIG. 1: Expectation values of averaged potential term for the 
harmonic oscillator in the conditions A, B, and C (see text). 
Dashed line denotes the potential expectation value of the 
ground state (w/4). 



(B) , iV sum = 200, A, = 3, and the fixed scale with £ = 0.2 

(C) . For the initial condition we fix Tj as Tj = il~/N sum 
and randomly generate qi within [—A q , A q ] as a hot start. 
We prepare N = 400 paths for the three conditions, re- 
spectively. Since temperature of the system 1/7" is much 
smaller than the excitation energy, the quantum fluctu- 
ations in this study will reflect the ground state of the 
harmonic oscillator. 

As a result of the numerical simulation, quantum paths 
approach to equilibrium at around the iteration number, 
i.e., number of steps 2-3., -ZVjt era tion ~ 10 3 for the condi- 
tion A, Alteration ~ 3 X 10 3 for B, and Alteration « 9 X 10 3 

for C. In Fig. [T] we show cooling behavior of the po- 
tential expectation value Vho = muj 2 q 2 /2 in averaged 

y 

form [A = J dr A(t) /T] for the harmonic oscillator by 
the Metropolis method. As one can see, at the satura- 
tion point the potential expectation value reproduces the 
ground-state value (= w/4) with more than about 90% 
accuracy in all conditions A, B, and C. Indeed, after it- 
eration Alteration = 10 4 the expectation value becomes 
(Fho) = 0.213 ± 0.04 for the condition A, 0.237 ± 0.004 
for B, and 0.241 ±0.004 for C. The potential expectation 
value approaches to the ground-state value as the scale 
constant £ gets small and hence structures of the quan- 
tum fluctuations in fine scale can be described, which is 
similar to the case of the discretized approach. We ex- 
pect that the accuracy of the potential expectation value 
will get higher if one uses even smaller £, although in 
which case large N SU m, large Alteration, and large simu- 
lation time are required. We note that the expectation 
values are saturated with respect to N sum , namely we ob- 
tain same expectation values within range of statistical 
errors with larger N sum in each conditions. 

In order to see the quantum fluctuations in detail, we 
show in Fig. [H upper) examples of the quantum fluctua- 
tions out of the N = 400 paths. From the figure, in all 
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FIG. 2: (Upper) Examples of quantum fluctuations for the 
harmonic oscillator in the conditions A, B, and C (see text). 
(Lower) Distributions of the coordinate for the harmonic os- 
cillator in logarithmic scale together with the squared wave 
function of the ground state denoted by dashed line. 



conditions the paths fluctuate from the origin q = to 
their maximal amplitude ~ ±1.5. The mean squared ra- 
dius of the fluctuation is (q 2 ) = 0.426±0.08 for the condi- 
tion A, 0.474± 0.008 for B, and 0.482± 0.008 for C, which 
arc close to the ground-state value (g 2 )cs = 0.5 within 
smaller than about 10% discrepancies. We note that in 
all conditions fluctuations in large scale are mainly com- 
posed of peak structures of width > 1 . Especially in the 
conditions B and C the fluctuation structures in large 
scale arc described by sum of several Gauss functions, 
although they have small fluctuations in fine scale. Then 
let us visualize degree of the quantum fluctuations. For 
this purpose we make a histogram for q with division 
of time range into sufficiently many parts in each path 
and then combine A" = 400 histograms to obtain the 
g-distribution p q . The result is shown in Fig. flower) 
together with the squared wave function of the ground 
state. As one can see, our ^-distributions in all condi- 
tions behave consistently with the squared wave func- 
tion. Especially it is interesting that behavior of the 
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quantum fluctuations to large q (~ ±3) is very similar 
to the squared wave function. 

In the examination of our approach for the harmonic 
oscillator, we have seen that our approach reproduces 
quantum behavior of the system with more than about 
90% accuracy by properly chosen scale constants, both 
in random and fixed scale cases. Especially, by using 
smaller scaling constant £, quantum behavior with higher 
accuracy is obtained. 

Next let us apply our approach to rclativistic field the- 
ories. It is important that our approach has possibilities 
to become a nonperturbative way to quantum field the- 
ories. To be specific, we firstly consider U(l) gauge field 
A fM (x) in four dimensions [x = (x, r), \x = 1, 2, 3, 4], 
and assume a periodic boundary condition with box size 
(X , X, X, T). In a similar manner to the nonrelativistic 
particles, the quantum transition amplitude of the field 
in Euclidean time can be expressed as, 



0.004 



Z = f T>Aexp(-S[A]), VA = JJdA^a;). 

./period Xjl 



(8) 



In the field path integrals, sole difference to the particle 
case is that the field is a function of four components 
of the coordinate x rather than time r only. Therefore, 
smooth quantum fields can be obtained by applying our 
approach (steps 1-4.) with an extension of r — > x. Then 
the field fluctuations are expressed as, 



A n( x ) = E Ai " exp 



(9) 



for each \i, where (x — Xi ) 2 means to take squared dis- 
tance between x and Xi^ in the periodic boundary con- 
dition. 

Now let us evaluate quantum fluctuations of the U(l) 
and SU(2) gauge fields in four dimensions, which actions 
are, 



S = / d 4 xC(x), 



(10) 



with Lagrangian densities, 
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(11) 



£u « = | E ~ dvA 

f- L 1 ly — 1 

£su(2)^E E (d^-d^ + g^^bcAlA , 

a—1 fl.U—1 \ b.C / 

(12) 

respectively. In this study we do not include gauge fixing 
terms nor the Faddeev-Popov ghosts in the Lagrangian 
densities. The SU(2) gauge field has self-interactions 
with coupling g, for which we take g = 3.5. In both U(l) 
and SU(2) gauge theories, we take the following simu- 
lation condition. Namely, we fix the coordinate a?i , tt) 



-0.004 



-0.008 



U(l) — B— 




sum — 

Linear 




Coulomb 






I LU [fi 


a 'A 
* 

.-' i 
.■ i 

(i t i i 





1.2 



0.8 ^ 

0.6 S 

0.4 



3 4 5 



0.2 





FIG. 3: Potential between fundamental representations for 
U(l) (left axis) and SU(2) (right axis) gauge fields. The U(l) 
potential is negatively shifted so as to fit the Coulomb poten- 
tial V(R) = —a/R with a « 1/137 denoted by dashed line. 
Dotted line denotes the linear potential V(R) = oR + b fitted 
to the SU(2) potential. 



to be on sites of four-dimensional lattice 7 3 x 14 divid- 
ing the box in same intervals with T = 2X . We take 
the fixed scale with £ = X/(Jy/n) = T/(Uy/n). Wc 
note that £ corresponds to the minimal scale of quantum 
field theories, as the lattice spacing a in the discretized 
framework. The cut-off for the fluctuation amplitude, 
Aa, is fixed as = 1.3 At first of the itera- 

tion Ai , j is randomly generated as a hot start. We 
prepare N = 50 paths for the U(l) and SU(2) gauge 
fields, respectively. During cooling by the Metropolis 
method, the action of the U(l) [SU(2)] gauge field con- 
verges at around Alteration ~ 3 x 10 5 (10 6 ). It is interest- 
ing that at the saturation point (£su(2)) ~ 0.49 £ -4 is 
smaller than 3 x (£tj(i)) ~ 3 x 0.20 due to the self- 
interactions in SU(2), where C is averaged Lagrangian 
density [Z = / d 4 xC(x)/TX 3 }. 

Quantum fluctuations of gauge fields provide a poten- 
tial between (infinitely heavy) fundamental representa- 
tion and its antiparticle, which can be evaluated through 
the Wilson loop of rectangle C = T x R defined as @ , 



W(T, R) =trVcxp 



Y,dx,A;(x)T a 



(13) 



where V means to take the ordered exponential with the 
group generator T a . We choose the gauge coupling for 
U(l) as g = 0.303 so that a = g 2 /( 47r ) ~ 1/137. From 
the expectation values of the Wilson loop (W), the po- 
tential is evaluated as, 



v ' t (W(T + t,R)Y 



(14) 



for sufficiently large T and small t. In this study, in 
order to have enough statistics, wc calculate average of 
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FIG. 4: Distribution of the U(l) field value A li (x) (open 
and closed circles) and distribution of the coefficient of the 
Gauss function for the U(l) field (open and closed squares), 
with different cut-offs for the fluctuation amplitude, Aa, in 
logarithmic scale. 



10 Wilson loops at random positions in four-dimensional 
time-space for each set (T, R) in each path and then com- 
bine the results of N = 50 paths. The results for U(l) 
and SU(2) gauge fields are shown in Fig. [3l As one can 
see, while the U(l) gauge field qualitatively reproduces 
the Coulomb force, the potential from the SU(2) gauge 
field shows confining linearity. A confining potential is 
the expected nonperturbative property of non-Abelian 
gauge theories inspired by the absence of free quarks in 
experiments [Tpj ] . and our approach implies that quanti- 
zation indeed generates confining field configurations in 
the SU (2) gauge theory. Fitting our SU(2) potential with 
aR + b we obtain the string tension a = 0.17 ± 0.01 £~ 2 , 
which brings a scale to the quantum SU(2) gauge theory. 

Here we should discuss the gauge symmetry in our 
approach. In the path integral formulations of gauge 
theories without gauge fixing terms nor Faddeev-Popov 
ghosts, the measure DA and the Lagrangian £ are re- 
spectively gauge invariant. This indicates that numeri- 
cal simulations in such a condition will have the gauge 
symmetry if one takes into account all the possible paths 
required by the measure, that is, the full gauge group 
manifold, which includes all of the gauge copies and es- 
pecially all of the Gribov regions in SU(2) [ll|]. In our 
approach, however, we take into account not all but most 
of the possible paths, as discussed above, and hence our 
approach does not contain the full gauge group manifold. 
Then we have two factors which leads to gauge symme- 
try breaking; one is that number of the Gauss functions 
N sum is finite and hence lack of the some paths may take 
place, and the other is the cut-off for fields, Aa, because 
we neglect regions out of the cut-off in the simulations. 
The first factor will become unimportant and negligible 
when one uses large value of N sum so that distribution of 
the Gauss functions is sufficiently dense. As for the sec- 
ond factor, on the other hand, regions out of the cut-off 



might contribute to the path integrals. In order to see 
the cut-off dependence for our results, we simulate quan- 
tum fluctuations of the U(l) and SU(2) gauge fields with 
the cut-off Aa = 2.6£ _1 , which is twice larger than that 
in our preceding simulation. We check that the cut-off 
dependence of the expectation values of the Lagrangian 
densities and potentials between fundamental represen- 
tations for both U(l) and SU(2) fields is negligible. Also 
we plot in Fig. [4] the distributions of the U(l) field value 
A^x) and of j4» , the coefficient of the Gauss function, 
with two cut-offs A^ = 1.3£ _1 and A^ = 2.6£ _1 . Here 
the distributions are evaluated in similar manners to the 
distribution of the coordinate for the harmonic oscillator 
discussed above. From the figure, both A fl (x) and A if ^ 
dominantly distribute around zero while the distributions 
for large values are suppressed at least exponentially be- 
cause of the weight exp(— S) in the simulations. Due to 
the exponential suppression, contributions of the gauge 
field configurations to the path integrals are effectively 
taken into account if one uses sufficiently large cut-off, 
which is the case in our simulations, and as a conse- 
quence the cut-off dependence cannot be seen for both 
A^x) and A ifi in Fig. [4] This is the result for the U(l) 
gauge field, but we also obtain a similar result for the 
SU(2) gauge field. In this sense, we consider an effective 
gauge group manifold than the full gauge group mani- 
fold, and in principle the effective manifold can be taken 
as close to the full manifold as possible by considering 
dense distribution of Gauss functions with large cut-off 
A a- Furthermore, this fact leads to a conjecture that 
in general quantum fluctuations of gauge fields appear 
dominantly within certain band and fluctuations out of 
the band is suppressed due to the weight exp(— S). 

In summary, we have developed an approach to eval- 
uation of the continuum path integrals, in which paths 
are described by sum of smooth functions with weight 
cxp(— S) by the Metropolis method. In this study we take 
an approximation that smooth fluctuations are described 
only by the Gauss function. The weighted smooth paths 
reproduce properties of the ground-state harmonic oscil- 
lator in one dimension with more than about 90% accu- 
racy by properly chosen width of the Gauss functions. 
We have found that quantum behavior with higher ac- 
curacy is obtained by using smaller width of the Gauss 
functions, with which finer structure of the quantum fluc- 
tuations can be described. We have also evaluated quan- 
tum fluctuations of fields and the Coulomb force and con- 
fining linear potential have been extracted at qualitative 
levels from the U(l) and SU(2) gauge fields in four di- 
mensions, respectively. We have found that distributions 
of large values of the gauge fields are suppressed at least 
exponentially in our approach, which implies that con- 
tributions of the gauge field configurations to the path 
integral are effectively taken into account by sufficiently 
large cut-off for the fluctuation amplitude. 
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Scientific Research (No. 22-3389). The author acknowl- 
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